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We present a calculation of the K — + tttt decay amplitudes from the K n matrix elements using leading 
order relations derived in chiral perturbation theory. Numerical simulations are carried out in quenched QCD 
with the domain- wall fermion action and the renormalization group improved gluon action. Our results show that 
the 7 = 2 amplitude is reasonably consistent with experiment whereas the 7 = amplitude is sizably smaller. 
Consequently the A7 = 1/2 enhancement is only half of the experimental value, and e' /e is negative. 



1. INTRODUCTION 

Quantitative understanding of the K tttt de- 
cay including the A7 = 1/2 rule and the value 
of e'/e has been a long-standing issue in lattice 
QCD lit]. Here we present a summary of our 
work H based on (i) the leading order chiral 
perturbation theory (ChPT) formula relating the 
K —^ TTTT decay amplitudes to the X — > tt matrix 
elements ^J, (ii) domain- wall fermion action to 
ensure chiral symmetry at finite lattice spacings 
in principle, and (iii) the renormalization group 
(RG) improved gauge action for suppression of 
chiral breaking effects for moderate fifth dimen- 
sional lattice size over those for the standard pla- 
quette gauge action H| . For a similar work using 
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Table 1 

Numbers of configurations generated in our cal- 
culation for each combination of and volume. 



ruf 


0.02 0.03 0.04 0.05 0.06 


16^ X 32 
243 X 32 


407 406 406 432 435 
432 200 200 200 200 



the plaquette gauge action we refer to Ref. 



2. NUMERICAL SIMULATION 



The numerical task at 
late the matrix elements 



hand 



Q 



to 



calcu- 
and 



1(0) 



where Q['\i = 1,---,10) are the 
local 4-quark operators of isospin / = 0, 2 that 
appear in the effective weak Hamiltonian. 

This calculation is carried out with the RG- 
improved gauge action at (3 = 2.6, which corre- 
sponds to the scale 1/a = 1.94 GeV determined 
from the string tension y/a — 440 MeV. For the 
domain-wall fermion action, we use the domain 
wall height M = 1.8 and the fifth dimensional 
size Ns = 16. We consider the case of degen- 
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Figure 1. 7r+ matrix element of Qe"'' for 

(a) 16'^ and (b) 24^. Solid and dashed curves show 
chiral fits as described in text. 

erate u-d and s quark with a common mass mj. 
Gauge configurations are generated independently 
for each value of m j . Our statistics are summa- 
rized in Table |^ for the two lattice volumes and 
m/ used in the calculation. 

3. CHIRAL PROPERTIES 

Our representative results for the X — > vr ma- 
trix elements are shown in Figs. ^ and ^j. In terms 
of the lattice pion mass squared, m^^, the — > tt 
matrix elements are expanded as 



Q 



K 



= flo + aim|f + a2{m\jY 



(1) 



Checking the chiral property with a quadratic fit 
with the fit parameters (09,01,02), we find the 
constant oq to be consistent with zero within sta- 
tistical errors for most of matrix elements. The 
only exception is Q^^^ as seen in Fig. How- 
ever, fits with the constraint oq = such as a 
cubic one using (01,02,04) or a form with the 
logarithm (01,02,03) yield acceptable x^/DOF. 
We therefore consider that the chiral properties 
required for viability of the ChPT relations be- 
tween the K — > TTTT and K n amplitudes hold 
well. 

An open problem is a small value of the coeffi- 
cient of the chiral logarithm term obtained in the 
fits, e.g., 03/01 = —0.58 GeV~^ for Qp\ which 
is only about 25% of the prediction of chiral per- 
turbation theory [§: -3/(87r2/^)- 
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Figure 2. Same as Fig. |l| but for Q^i\ 

4. PHYSICAL RESULTS 

We construct the K tttt decay amplitudes 
from the K ^ tt matrix elements as described 
in Ref . ||] . Renormalization of the lattice matrix 
elements is made with tadpole-improved pertur- 
bation theory at one-loop level at the match- 
ing point q* = 1/a. The renormalization group 
running from fi — 1/ a to rUc (=1.3 GeV) is car- 
ried out for the Nf = 3 theory with — 325 
MeV Combining these results with the Wil- 
son coefficients lead to the physical results for the 
K ^ TTTT decay amplitudes Aj {I = 0, 2). 

In Fig. ^ we present ReA2, KeAo and lo"^ = 
KeAo/KeA2 as a function of m\.j. We observe 
that the data from the volumes 16"^ (open sym- 
bols) and 24'^ (filled symbols) are consistent. 
Lines drawn are the fit curves for the quadratic 
(solid) function in or a form including chi- 
ral logarithm (dashed). The / = amplitude 
reaches only half of the experimental value in the 
chiral limit, and hence the A J = 1/2 rule, i.e., 
uj~^ « 22, is reproduced only by about 50% in 
the chiral limit. 

In Fig. I we plot p(3/2) p(i/2)^ which are 
related to e'/e as H as 



(2) 



Here we use our results only for the factor Im^o,2- 
substituting the experimental value for the real 
part of the amplitudes. Our results for e' /e is 
negative since p(3/2) ^ p(i/2) ^ 

5. PENGUIN OPERATORS 

Recently, Golterman and Pallante |^ pointed 
out a necessity of new operators to represent the 
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Figure 3. Real parts of the decay amplitude KeAj [10~*GeV] and the ratio ui^^ = Re^o/R-e^2 as a 
function of m\j [GeV^]. Open (filled) symbols are from the volume IB'^ (24^), and solid (dashed) curve 
show quadratic (chiral logarithm) fit result. 
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Figure 4. Imaginary parts of the decay amplitude P^^/^\ p(i/2) g^y^^ e'/e. Meaning of symbols is same 
as in Fig. ^. 



penguin operators in (partially) quenched chiral 
perturbation theory, and gave a simple prescrip- 
tion to extract the physical contribution of Qq^ 
to the if — > TT and K ^ matrix elements in this 
framework. Our preliminary analyses show that 
this modification, if employed, increases 

p(l/2) by 

about 10-50% depending on m/, but this effect 
is too small to affect the negative value of e'/e we 
have obtained. 
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